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Confexf. As increasingly more spectroscopic data are being delivered by medium- and high-resolving power multi-object spectro- 
graphs, more automatic stellar parameter determination softwares are being developed. The quality of the spectra collected also allows 
the determination of elemental abundances. 

Aims. SPADES is an automated software for determining: the radial velocity (Vr), the effective temperature (r cff ), the surface gravity 
(logg ), the metallicity ([Fe/H] ), and most importantly, the individual abundances. In this first version it is targeted on the analysis 
of mid-F-G dwarfs, but is meant to evolve to analyze any type of single stars. 
Methods. SPADES relies on a line-by-line modeling to determine the stellar parameters. 

Results. The internal systematic and random errors of SPADES were assessed by Monte Carlo method simulations with synthetic 
spectra and the external systematic errors by analysing real ground-based observed spectra. For example, by simulating the Giraffe 
setups HR13 and HR14B with synthetic spectra for a dwarf with T cf j = 5800 K, logg = 4.5, [Fe/H] = 0.0 dex and with a signal-to- 
noise ratio (snr) of 100, the stellar parameters are recovered with no significant bias and with l-cr precisions of 8 K for T e g, 0.05 for 
log g , 0.009 for [Fe/H] , 0.003 for [Ti/Fe] and 0.01 for [Ni/Fe] . 
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1. Introduction 

The present and future large spectroscopic surveys are going 
to significantly increase the number of spectroscopic data to 
be analyzed. A few e xamples are RAVE with some 250 00 
stars observed so far dBoeche et al.ll201 ll ISiebert et ai1l201 ll) . 
the Gaia-ESO Survey with ab out 10 s stars to be observed, Gaia 
with about 200 million stars dKatz et al.ll2004T: lKatzll20"09l) . To 
analyze these quantities of data, automatic spectra analysis 
softwares are needed. 

The development of automated spectra analysis software can 
be traced back to the early 1990s. Since then, many programs 
were developed based o n a wide variety of methods: e.g. ortho 
onal vector projection (ICavrel et all 1 1 99 U iPerrin et al.l Tl995 
minimum distance method, e.g. TGMET ( Katz et al] 1998; 
ISoubiran et al.ll2003l). and ETOILE dKatzll200ll). artifici al neural 
networks (Bailer- Jones 2000; Allende Prieto et al. 2000), feature 
fitting softwares such as Abbo dBonifacio & Caffau 2003) , pro- 
jection vector softwares su ch as MATISSE dRec io-Blanc o et alJ 
2006; Biiaoui etaTJ 120081). principal c omponent analysis soft- 
wares such as MAy dJofre etai]|2010l). and forward modeling 
algorithm ILIUM (Bailer- J ones! 120101). One can also cite 
DAOSPEC dStetson & Pancinol I2008L l2010h and FITLINE 
dFrancois et al.ll2003l) . which automatically measure equivalent 
widths. DAOSPEC and FITLINE are not automated per-say 
parameter determination softwares, but coupled with programs 
that derive atmospheric parameters from equivalent widths, they 
are used by many astronomers to automatically parameterize 
stars. 

The work presented here describes the development of a 
new automatic stellar spectra analysis software. This is based 



on line-by-line analysis. The reference lines are modeled. One 
of the motivations to develop SPADES (Stellar PArameters 
DEtermination Software) is the analysis of a set of 200 Giraffe 
HR13 and HR14 spectra, mainly F-G dwarfs, that were collected 
in the context of a study of the structure of the Galactic disk. 
Therefore, the first version of SPADES is illustrated and tested 
with simulated Giraffe-like spectra of mid-F-G dwarfs. SPADES 
is coded with Java. 

SPADES is described in Sect. [2] It relies on auxiliary data: 
a line list and a grid of reference spectra. They are presented in 
Sect. [3] The performances of SPADES were evaluated by Monte 
Carlo method simulations against synthetic spectra and by com- 
parison to observed spectra of the Sun, v And, (3 Vir, jj Her, and 
<t Dra (see Sect. |4j. 



2. SPADES 

2.1. General concept 

SPADES automatically determine: the radial velocity (Vr), the 
effective temperature (r e ff), the surface gravity (logg ), the 
metallicity ([Fe/H] ), and individual abundances ([X/Fe] for 
the element X) for single stars. This first version of SPADES is 
focused on mid-F-G dwarfs and the micro-turbulence is set to 
1 kms -1 . 

Many of existing automatic spectra analysis softwares uses 
global methods, meani ng they proc ess a ll pixels of the stud- 
ied spectrum dKatz et al.l 1998UKat3l200lHSoubiran et alJl2003b 
iRecio-Blanco et al.ll2006l) . With SPADES we aimed to explore 
another possibility, which is the line-by-line analysis. Each pa- 
rameter is determined using a pre-defined set of lines. The merits 
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of this technique are, first, the possibility to select lines that are 
particularly sensitive to the stellar parameter to be determined, 
and second, this allows one to exclude lines that are not correctly 
modeled. A drawback is that only a fraction of the total available 
information is used. 

Often, whether automated or manual, methods using a line- 
by-line analysis rely on equivalent width measurement and use a 
curve-of-growth analysis. We chose to rely on synthetic spectra 
modeling and profile fitting. The motivation is twofold. First, this 
takes care of small blend problems. Assuming that the blend is 
correctly modeled, its contribution to the observation and to the 
modeled profile will (partly) cancel out. Second, some parame- 
ters can be derived by line profile fitting, e.g. r e fffrom the wings 
of the Balmer lines. A drawback of this method is the need for 
a very large library of synthetic spectra. Reading and process- 
ing many synthetic spectra is computationally demanding and 
therefore the analysis is slower. 

In the line-by-line analysis the accuracy and precision of the 
results is very sensitive to the determination of the continuum or 
pseudo-continuum. Yet, this is often difficult particularly in cool 
and/or super metal-rich stars. SPADES determines "local" con- 
tinuum from the comparison of the reference (synthetic) spectra 
and the studied spectrum (see sect. 12.5.31 ). 

As opposed to mean metallicities, SPADES determines iron 
and individual elements abundances. To derive of the individual 
abundances SPADES calculates on-the-fly atmospheric models 
and synthetic spectra for the values of T e ff, \ogg , [Fe/H] it had 
determined (see sect. 12.61 ). 

SPADES can simultaneously analyze various spectral do- 
mains (possibly with different resolving powers) of the same tar- 
get, e.g. different Giraffe setups, or different orders of an echelle 
spectrograph, or even spectral domains collected with different 
spectrographs. 

2.2. Processing overview 

A SPADES analysis proceeds in three steps: it first determines 
the radial velocity, then the atmospheric parameters, and finally 
it determines the individual abundances. 

The radial velocity is determined by cross-correlation with a 
synthetic template (see sect. I2.41 i. 

The three atmospheric parameters (effective temperature, 
surface gravity, iron over hydrogen ratio) are determined inde- 
pendently (see sect. 12.51 ). Currently, each parameter is associated 
to one diagnostics. SPADES was conceived to aggregate more 
diagnostics along its life span and more will be added in the fu- 
ture Current diagnostics are 

- effective temperature: Balmer line wings fitting (suited for 
r e fjderivation in F-G dwarfs), 

- surface gravity: ionization equilibrium, 

- iron abundance: neutral iron line profile fitting. 



The individual abundances 
element by line profile fitting. 



are determined element per 



The reader will notice that some of these diagnostics are sim- 
ilar to those used in the so-called "detailed analysis". Two differ- 
ences should be emphasized: (i) SPADES relies on profile fitting 
and modeling rather than equivalent width measurement, and (ii) 
the SPADES precision is not limited by the step in the grid of at- 
mospheric models and spectra (see Sect. l2.5.TT l. 



2.3. Pre-processing 

Each of the three processing steps includes pre-processing 
tasks. Pre-processing a spectrum includes shifting to rest frame, 
excluding cosmics, excluding telluric lines, convolving the 
synthetic spectrum to the resolving power of the ground-based 
observed spectrum, normalizing to continuum and excluding 
possible unexposed edges of the spectrum. Every functionality 
remains an option: it can be switched on and off in the configu- 
ration file by the user. 

Shifting to rest frame is performed by using either a radial 
velocity provided by the user or that determined by SPADES 
(see sect. 12.41 ). 

The exclusion of cosmics is performed by detecting the 
fluxes that lie significantly above the continuum and replacing 
their value by interpolation with the neighboring pixels. 

The wavelengths of the telluric lines were pre-tabulated. 
These lines are not cut. Their pixels are simply flagged as not 
valid and are excluded from all subsequent analysis. 

The spectral resolution of the synthetic spectra should be 
matched to the resolution of the ground-based observed spec- 
trum. To do so, the reference spectra are convolved using a 
Gaussian profile (valid for slow rotators). 



Raw spectrum 



Wavelength (A) 



Final spectrum 




6500 6550 
Wavelength (A) 



Fig. 1. Observed Giraffe/VLT (HR14B setup) Sun spectrum be- 
fore normalization, with the polynomial fit shown as a thick blue 
line (top) and after normalization (bottom) 



The synthetic spectra are already normalized to continuum. 
The studied spectrum is iteratively fitted with a polynomial func- 
tion, using a sigma-clipping technique to exclude the lines and 
cosmic rays from the fit. The degree of the polynomial is defined 
by the user, depending on the large-scale shape of the spectrum. 
The studied spectrum is normalized to 1 at the level of the con- 
tinuum by division by the polynomial. Figure Q]presents an ob- 
served Giraffe/VLT (HR14 setup) Sun spectrum before normal- 
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ization, with the polynomial shown as a thick blue line (top) and 
after normalization (bottom). 

Edge exclusion is necessary if the ground-based observed 
spectrum suffers from underexposed or unexposed pixels. 
SPADES can automatically detect and flag as not valid the un- 
derexposed/unexposed pixels. 

2.4. First step: Radial velocity 

The radial velocity determination is based on a cross-correlation 
of the studied spectrum and a synthetic template spectrum. 

The atmospheric parameters of the template spectrum are 
chosen by the user. At this stage in the processing, SPADES 
has not yet provided estimates for the atmospheric parameters 
of the studied source. If no external information is available 
to guide the choice of the template, it is likely that there will 
be a mismatch between the ground-based observed spectrum 
and the selected template. Monte Carlo method simulations 
were performed to assess the impact of selecting a template 
400 K and 1 000 K warmer than the studied spectrum. The tests 
were conducted with 200 snr=50 synthetic studied spectra of 
solar parameters. In both cases (400 and 1 000 K offsets), the 
mismatch leads to a bias of about 0.2 kms" 1 and no significant 
change in the precision of the estimated radial velocities (see 
table [TJ. A bias of 0.2 kms -1 represents about l/60' /! of the 
full width at half maximum (FWHM) and has a negligible 
impact on the following steps of the analysis. If, a posteriori, 
the parameters estimated by SPADES and those of the selected 
template are significantly different, the radial velocity can be 
re-estimated with a better matching template. 



Correlation function 
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Fig. 2. Correlation function for a synthetic studied spectrum 
of solar parameters (T e ff = 5777 K, logg — 4.44 and [Fe/H] 
— 0.0 dex) and S/N =100 (Poisson noise) and a noiseless syn- 
thetic template with the same atmospheric parameters 



The studied and the template spectrum are pre-processed. 
The reference spectrum is then shifted step by step in radial ve- 
locity over an interval defined by the user. In the tests, we used 
an interval of -500 to 500 km/s. At each step, the synthetic spec- 
trum is re-sampled to the same sampling as the studied spectrum. 
At each step, the correlation coefficient of the two spectra is cal- 
culated: 



This allows one to estimate the radial velocity with a precision 
better than the step used to shift the template. 

The error on the radial v elocity determination, cr Vr is esti- 
mated with the Zucker (2003) formula: 



CC(v) = - J] - ° V K S > > X lRU) ~ A ' >> 



!=0 



stdv(S) x stdv(R) 



(1) 



Vr 



NX 



C"(m) C(m) 1 



C(m) 1 - C{ m ) 



(3) 



where v is the velocity shift applied to the template spectrum, 
N is the number of valid pixels in the studied spectrum, S (resp. 
R) the studied (resp. reference) spectrum counts, < S > (resp. 
< R >) the mean studied spectrum (resp. reference) counts, and 
stdv the standard deviation: 



1 N 

stdv(X) = (^2 (X(0_ < X >)2 ) 



1/2 



(2) 



The correlation coefficients as a function of the radial veloc- 
ity shifts form the cross-correlation function. Figure|2]shows the 
cross-correlation function obtained for a synthetic studied spec- 
trum of solar parameters (T e g = 5777 K, log g — 4.44 and [Fe/H] 
= 0.0 dex) and S/N = 100 (Poisson noise) and a noiseless syn- 
thetic template with the same atmospheric parameters. 

The radial velocity of the observed source corresponds to 
the maximum of the correlation function. The maximum is 
derived in two steps. First, the highest correlation coefficient is 
found and then a second-degree polynomial is adjusted to the 
correlation function in the vicinity of the highest coefficient. 



where C(m) is the value of the correlation function at its maxi- 
mum, C"(s) the value of the second derivative of the correlation 
function at the maximum of the correlation function and N the 
number of valicQ pixels in the ground-based observed spectrum. 

The performances of the radial velocity estimation of the 
radial velocity were assessed with a Monte Carlo method. For 
three S/N: 30, 50 and 100, Giraffe HR13 and HR14B spectra 
were simulated (200 per SNR) by injecting a synthetic spec- 
trum with Poisson noise. The synthetic spectrum was a solar-like 
spectrum in the rest frame. The same spectrum, but noiseless, 
was used as template. Table Q] presents the median and disper- 
sion of the distributions of radial velocities derived by SPADES 
for each SNR. Because the true radial velocity of the studied 
spectra is kms~' , the median of a distribution is an estimate 
of the systematic error and the dispersion an estimate of the ran- 
dom error. 



1 Some non-valid pixels affected by, e.g. telluric lines, could be ex- 
cluded from the derivation of the correlation function (see Sect. 12.31 , 
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Table 1. Median and dispersion (in km s 1 ) of the radial velocity 
distribution derived by SPADES for 200 studied spectra for SNR 
30, 50 and 100. 



SNR 


Vr median 


a Vr 




(km s" 1 ) 


(kms- 1 ) 


100 


-0.005 


0.03 


50 


-0.006 


0.06 


30 


0.002 


0.11 



2.5. Second step: Stellar parameter determination 
2.5.1 . Overview 

SPADES determines the effective temperature (7W), the sur- 
face gravity (\ogg ), and the iron abundance ([Fe/H] ). In this 
first version targeted on mid-F-G dwarfs, we assumed a micro- 
turbulence of 1 kmr 1 for all our stars, including the tests with 
the observed spectra. 

The determination of the atmospheric parameters is iterative. 
Initial values for the three parameters are provided either by 
the user or by SPADES (see Sect. 12.5.2b . Each parameter is 
determined independently of the two others, leaving these two at 
their initial values. At the beginning of a new iteration, the val- 
ues of the three parameters are updated, therefore providing new 
"initial" values. This is repeated until the convergence condition 
is met: i.e. absolute differences between the parameters obtained 
at iteration n and n — 1 all simultaneously lower than 10 K for 
T e g, 0. 1 for log g , and 0.025 for [Fe/H] . The iterative approach 
was chosen over a direct approach (e.g. scanning the whole grid 
of 2000 reference spectra over all parameters at once), because 
it was computationally faster. Yet, the iterative approach could 
potentially be more complex in terms of convergence. The 
convergence performances of SPADES are assessed in Sect. 14.31 

At any iteration, the estimation of each parameter follows the 
same logic: 

SPADES first calculates a ID grid of synthetic spectra, varying 
only along the parameter that is estimated. This lD-subgrid is 
calculated from the large grid of synthetic spectra described in 
Sect. 13.11 The grid is linearly interpolated to the fixed values of 
the other parameters. As an example, at the start of iteration n, 
the starting values for the atmospheric parameters derived from 
iteration n-1 are r eff = 5743K, logg = 4.42, and [Fe/H] = 
-0.02. The first parameter to be determined in iteration n is T e g. 
A lD-subgrid is linearly interpolated for r e (j ranging from 5200 
to 6200 K by step of 200 K and logg = 4.42 and [Fe/H] = 
-0.02 dex. 

Each synthetic spectrum is then compared (not globally, but 
the selected lines) to the ground-based observed spectrum ac- 
cording to the diagnostic appropriate for the parameters to be 
determined. As an example for T e ff, the Ha profile of the ground- 
based observed spectrum will be compared to that of each syn- 
thetic spectrum. 

The result of the comparison is a similarity indicator for each 
node in the lD-subgrid. This series of similarity indicators as a 
function of the parameter to be determined form a function that 
is analyzed (minimized or nulled, depending on the function) 
for the best similarity value, leading to the best estimate of 
the parameter to be determined. Interpolation or polynomial fit 
allows one to determine the parameter with a precision that is 
not limited by the step of the lD-subgrid. 



SPADES is iteratively searching for the best matching syn- 
thetic spectrum in a 3D space (r e ff, logg and [Fe/H] ). In this 
first version of SPADES, the search is split in three ID minimiza- 
tions, rather than using a multi-dimensional search method. The 
motivation is that this facilitates visualizing and interpreting the 
residuals and similarity functions (See e.g. Fig. [5]) that depend 
on a single parameter. In future versions, it would be interesting 
to experiment with multi-dimensional minimization algorithms 
such as Gauss-Newton, gradient, or Levenberg-Marquardt meth- 
ods. 

2.5.2. Starting parameters 

SPADES proceeds iteratively to derive the atmospheric param- 
eters. It therefore needs starting values for the effective temper- 
ature, surface gravity, and iron over hydrogen ratio prior to the 
first iteration. Those can be deduced from photometric or bibli- 
ographic information if available. 

If no external information is available, SPADES offers the 
possibility to perform a first parameterizati on of the studied 
source, using the TGMET/ETO ILE method dKatz etalj[l998t 
lKatzll200lt ISoubiran et aT]l2003h . The results of this first anal- 
ysis are used as starting values for the iterative parameterization 
with SPADES. 

The TGMET/ETOILE method is direct (not iterative). It 
compares a studied spectrum to a library of reference spectra, 
looking for those most similar in the least-squares-sense. The 
parameters deduced for the studied spectrum could either be the 
parameters of the most similar spectrum or a combination of the 
parameters of the most similar spectra. Because SPADES is only 
looking for starting values, it simply adopts the parameters of the 
most similar spectrum. 

The performances of the SPADES implementation of the 
TGMET/ETOILE method (hereafter referred to as SPADES- 
TE)) were assessed with a Monte Carlo method method. Starting 
from a synthetic spectrum with T e g = 5800 K, logg = 4.5, 
[Fe/H] = -1 dex (which is a node of the grid of synthetic spec- 
tra), three series of 200 studied spectra each were simulated for 
S/N=100, 50 and 30, respectively, and the resolving power and 
wavelength range of the Giraffe HR13 and HR14B setups. At 
S/N=50 and 100, the correct node in the grid is recovered in 
100% of the cases. At S/N=30, SPADES-TE returns most of the 
time a solution biased by one grid step in r e ff(i.e. +200 K) and 
two grid steps in [Fe/H] (i.e. +0.4 dex). 

SPADES-TE was also tested on a Sun spectrum observed 
with Giraffe HR13 and HR14B setups in twilight condi- 
tions. SPADES-TE found T eS = 5800 K, logg = 3.5 and 
[Fe/H] = -0.2 dex, to be compared to the canonical parame- 
ters for the Sun T eff = 5777 K, log g = 4.44, and [Fe/H] = dex. 

We recall that the version of SPADES presented here is lim- 
ited by the GRID of synthetic spectra and by the list of reference 
lines to the mid-F-G dwarf stars. Future versions of the program 
will rely on larger grid(s) and on several fine-tuned line lists, 
which will extend the scope of SPADES. 

2.5.3. Line and continuum identification and vertical 
adjustment 

For each spectral line, the pixels constituting the line are 
dynamically identified using the synthetic spectra. Using the 
rest wavelength of the line, SPADES starts from its center and 
agregates pixels in each wing as long as the value of the signal 
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6290 6295 6300 
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Fig. 3. Two examples: (i) pixel identification in the reference line 
(green crosses), (ii) pseudo-continuum selection (blue circles), 
and (iii) vertical adjustment of the synthetic reference spectrum 
(red line) to the studied spectrum (black line). In the top plot, 
the pseudo-continuum area (blue circles) were selected automat- 
ically by SPADES and were used for the vertical adjustment. In 
the bottom plot, no valid pseudo-continuum area was found. The 
pixels in the line (green crosses) were used for the vertical ad- 
justment. 

increases. If less than five pixels are associated to the line, it is 
considered too weak and rejected from the subsequent analysis. 



This parameter is adjustable by the user. Figure [3] shows two 
examples of sets of pixels associated to a line, i.e. green crosses 
on the figures. 

SPADES relies on profiles fitting. This requires adjusting the 
observed and synthetic spectra on the same vertical scale. It is 
possible, in the pre-processing, to normalize the overall contin- 
uum to 1. This is a global adjustment on the whole spectrum. 
This is not enough for the comparison of two lines. Therefore, 
SPADES systematically adjusts locally, line by line, the vertical 
level of the synthetic line to the one of the observed line. 

SPADES first identifies flat areas (pseudo-continuum) in the 
synthetic spectrum in the vicinity of the studied line and selects 
one on each side of the line. It then checks that these two areas 
also correspond to pseudo-continuum in the studied spectrum. If 
the check is positive, it uses these areas to vertically scale the 
synthetic spectrum on the studied one, by linear regression, i.e. 
there are two degrees of freedom, vertical scaling, and slope. 

If there are no flat areas in the synthetic spectrum or the 
studied spectrum contains (a) line(s) in the selected pseudo- 
continuum, another method is used for the vertical adjustement. 
The method depends on the diagnostic. For all diagnostics re- 
lying on the quadratic sum of the residuals (e.g. the iron abun- 
dance, Sect. 12.5.6b . the pixels in the lines are directly used for 
the vertical adjustment (allowing here for a single degree of free- 
dom, the vertical scale). For the diagnostic that relies on the the 
sum of residuals (i.e. the surface gravity, Sect. 12.5.51 . SPADES 
uses pre-tabulated pseudo-continuum. 

An exception to this scheme is the H a line used to derive of 
the effective temperature, for which the pseudo-continuum and 
line area are pre-tabulated. 

Figure [3] shows the local vertical adjustement of a synthetic 
spectrum (red) to a studied spectrum (black). In the top plot, 
the pseudo-continuum area (blue circles) were selected automat- 
ically by SPADES and were used for the vertical adjustment of 
the synthetic reference spectrum (red line) to the studied spec- 
trum (black line). In the bottom plot, no valid pseudo-continuum 
area was found. The pixels in the line (green crosses) were used 
for the vertical adjustment. 

2.5.4. Effective temperature 

In this article, the stars considered are mid-F-G dwarfs, the 
wings of t he Balmer lines are very sensitive to the effective tem- 
perature dvan't Veer-Menneret & Megessierl Il996t [Fuhrmannl 
119981 12004 12008t iBarklem et alj|2002t ICavrel et all 1201 lb . In 
this study, SPADES used the wings of H a present in Giraffe 
setup HR14B to derive the effective temperature. The core of 
H a was excluded because it is less sensitive to r e ffVariation and 
poorly modeled, specially using ID LTE models, which is the 
case here. 

As presented in section 12.5.11 a ID grid of synthetic spectra 
was calculated, differing only by the effective temperature. After 
pre-processing and vertical adjustment, each synthetic spectrum 
was compared to the studied spectrum. For each synthetic spec- 
trum, the degree of similarity of the synthetic H a wings with the 
studied one, is calculated as the quadratic sum of the residuals 
between the pixels in the two lines: 

N 

Sim Ha = J](S(i)-R(i)) z , (4) 

1=0 
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Teff function 



log(g) determination Function 




6540 6560 6580 

Wavelengths(A) 




6550 6555 6560 

Wavelengths(A) 

Fig. 4. Comparison of the H a wings in the studied spectrum 
(black line) with the ID grid of synthetic spectra (red lines) 
ranging here from 5000 (narrowest wings) to 6600 K (broad- 
est wing). The pseudo-continuum area is labeled with blue dots 
and thehave been valid pixels in the H a wings are identified with 
green crosses. The bottom plot is a zoom on the blue wing. 



where N is the number of valid pixels in the H a line, S(i) (resp. 
R(i) the counts of the studied spectrum (resp. reference spec- 
trum) for pixel i. 

Figure [4] shows the comparison of the H a wings in the 
studied spectrum (black line) with the ID grid of synthetic 
spectra (red lines) ranging here from 5000 (narrowest wings) to 
6600 K (broadest wing). The pseudo-continuum area is labeled 
with blue dots and the valid pixels in the H a wings are identified 
with green crosses. The bottom plot is a zoom on the blue wing. 

The degrees of similarity versus the effective temperature 
define a function that is minimum for the effective temperature 
of the studied source. The minimum is found by fitting a 
second-degree polynomial around the lowest similarity value. 
The top left plots in Figure [5] and [6] give examples of similarity 
functions as a function of effective temperature, for a synthetic 
spectrum and for the ground-based observed spectrum of 
v And respectively. Each synthetic spectrum in the ID subgrid 
corresponds to a red dot. The second-degree polynomial fitted 
to the similarity function to minimize it is represented as a blue 




5000 5500 6000 6500 
Tefl (K) 



[Fe/H] determination Function 



2.0 2.5 3.0 3.5 4.0 4.5 5.0 

l°9(8) 



Ni abundance determination Function 





-2.0 -1.5 -1.0 -0.5 0.0 
[Fe/H] 



-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 
[Nl/Fe] 



Fig. 5. Examples, for a synthetic studied spectrum, of similarity 
functions as a function of effective temperature (top left), sur- 
face gravity (top right), iron over hydrogen ratio (bottom left) 
and nickel over iron ratio (bottom right). The polynomial (or lin- 
ear interpolation for the surface gravity) used to minimize the 
functions are represented as blue lines. 



line. 



The similarity functions for the synthetic spectra and v And 
show similar degrees of regularity and smoothness. The first rea- 
son is that the spectrum of v And has a good S/N ratio. The sec- 
ond reason is that in each point (red dots) the similarity func- 
tions derives from the comparison of the same noisy studied 
spectrum with different noiseless synthetic spectra. The noises 
on the points of a similarity function are correlated. The noise in 
the studied spectrum and the systematic differences between the 
studied and synthetic spectra propagate to the similarity func- 
tions, not so much as fluctuation and irregularities in the func- 
tions, but as a vertical and horizontal shift and as a flattening of 
the curvature of the functions. 



2.5.5. Surface gravity by ionization equilibrium 

The ionization equilibrium method is based on the fact that the 
abundances determined from lines of the same element, but with 
different ionization stages, should be the same. Building on this 
idea, SPADES simultaneously models Fe lines in two ionization 
stages: Fel and Fell. 

Similarly as for T e ^, a ID grid of synthetic spectra, differ- 
ing only by the surface gravity (this time) was calculated. For 
each a "degree of ionization equilibrium" (hereafter refered to 
as "gravity similarity") is calculated as 



Sim„ 



j Nfel y N Fell 

= Y Res, V, 



Rest, 



(5) 
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Fig. 6. v And similarity functions as a function effective temper- 
ature (top left), surface gravity (top right), iron over hydrogen 
ratio (bottom left), and nickel over iron ratio (bottom right). The 
polynomial (or linear interpolation for the surface gravity) used 
to minimize the functions are represented as blue lines. 



where Nf c i (resp. NFeii) are the numbers of Fel (resp. Fell) lines 
used and R; the count residuals for the line 1: 



fo^ = ]T(S(0-tf(0), 

i=0 



(6) 



where N is the number of valid pixels in the line, S (f) (resp. R(i) 
the count of the studied spectrum (resp. reference spectrum) for 
pixel i. 

Figure [7] shows an Fel line (top, green crosses) and an Fell 
line (bottom, green crosses), used by SPADES to determine of 
log g . The studied spectrum (log g - 4.5) is denoted in black and 
the red lines are the synthetic reference spectra of the ID subgrid 
ranging from logg = 2.0 to \ogg = 5.0 by steps of 0.5. The neu- 
tral iron line is alsmost insensitive to gravity, while the singly 
ionized iron line becomes stronger for decreasing gravities. 

The gravity similarity, as defined above, should be null for 
the surface gravity of the studied source. The zero of the grav- 
ity similarity function is found by interpolation. The top right 
plots in Figure|5]and figure[6]give an example of similarity func- 
tions as a function of surface gravity. Each synthetic spectrum 
in the ID subgrid corresponds to a red dot. The linear interpo- 
lation used to find the gravity corresponding to the zero gravity 
similarity is represented as a blue line. 

2.5.6. iron over hydrogen ratio 

The [Fe/H] determination is based on a global Fel line profile 
fitting. For the other parameters, a ID grid of synthetic spectra, 
differing only by the iron abundance, was calculated. For each 




1 

6320 


6321 


6322 6323 
Wavelengths(A) 

Fell line 


6324 


I 

6325 






















6245 


6246 
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6249 


6250 



Wavelengths(A) 

Fig. 7. Fel line (top, green crosses) and an Fell line (bottom, 
green crosses), used by SPADES to determine of logg . The 
studied spectrum (logg = 4.5) is given in black and the red lines 
are the synthetic reference spectra of the ID subgrid ranging 
from log g = 2.0 to 5.0 by step of 0.5. 



synthetic spectrum, the global similarity between all synthetic 
neutral iron lines and all observed neutral iron lines is derived as 



Sinif e 



Y.Rej, 



1=0 



(7) 



where NFei is the number of neutral iron lines in the list of ref- 
erence lines, Resjis the square of the count residuals for line 1, 
given by formula|6] 

The iron content of the source corresponds to the minimum 
similarity degree. It is found by adjusting a third-degree polyno- 
mial through the "iron similarity function". A third-degree poly- 
nomial is preferred over a second-degree polynomial, because 
the "iron similarity function" is usually asymmetric: steeper on 
the metal-rich side than on the metal-poor side. The bottom left 
plots in Figure|5]and figure[6]give an example of similarity func- 
tions as a function of the iron over hydrogen ratio. Each synthetic 
spectrum in the ID subgrid corresponds to a red dot. The third- 
degree polynomial fitted to the similarity function to minimize it 
is represented as a blue line. 

2.6. Third step: Individual abundances 

Once the iterative determination of the three atmospheric param- 
eters has converged, the individual abundances are determined. 
The general principle is very similar to the determination of 
the iron content, except that the determination is direct and not 
within an iterative loop with several other parameters. 

Each element was determined individually, using a different 
pre-defined list of lines. The difference with the [Fe/H] deter- 
mination method is that for each element, the ID grid was cal- 
culated on-the-fly for the T^, \ogg , and [Fe/H] determined by 
SPADES and varying only by the abundance of the element to 
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be determined. SPADES interpolates on-the-fly the input files 
for the Kurucz programs (e.g. the opacity distribution functions 
- PDF (fCastelli & Kuruczl 12004 ) and calls Kurucz ATLAS9 
dKuruczll2005l) to calcula te the atmospheric models and Kurucz 
SYNTHE dKuruczll2005l) to compute the synthetic spectra. The 
GNU -Linux ported versions of ATLAS9 and SYNTHE were 
used (ISbordone et al.ll2004ISbordonell2005h . 

Element per element, for each synthetic spectrum the global 
similarity between all synthetic lines and all observed lines of 
the element X is derived as 

Sim x = ^ j Res 1 l , (8) 

z=o 

where Nx is the number of lines of the element X in the list of 
reference lines and Resjis the square of the count residuals for 
line 1, given by formula[6] 

The abundance of the element X corresponds to the minimal 
similarity degree. It is found by adjusting a third-degree poly- 
nomial through the element X similarity function. The bottom 
right plots in Figure[5]and figure[6]give an example of similarity 
functions as a function of nickel over iron ratio. Each synthetic 
spectrum in the ID subgrid corresponds to a red dot. The third- 
degree polynomial fitted to the similarity function to minimize it 
is represented as a blue line. 

3. Auxiliary data 

3. 1 . Reference grid of models and synthetic spectra 

The unidimensional reference grids used to derive the stellar pa- 
rameters are calculated by SPADES by interpolation using a pre- 
calculated grid. Models and synthetic spectra were calculated us- 
ing the Kurucz program s ATLAS9 and SYNTHE (lKuruczl2 005) 
ported to GNU-Linux (ISbordone etalJl2004t ISbordondl2005l) . 
About 2000 models and spectra were generated over the follow- 
ing range of parameters: 

- Effective temperature: r e fffrom 4800 to 6800 K with a step 
of 200 K. 

- Surface gravity: logg from 2.0 to 5.5 with a step of 0.5. 

- iron abundance: [Fe/H] from -3.0 to 1.0 dex with a step of 
0.20 dex. 

- Microturbulence: £ = 1 kms -1 

The models were calculated with the overshooting option 
switched off and l/H p = 1.25. The convergence of each 
model was checked. The spectra were calculated for the spec- 
tral domain ranging from 600 to 680 nm, corresponding to the 
Giraffe/VLT HR13 and HR14 setups and a resolving power 
R = 300 000 (SPADES convolves and re-samples the synthetic 
spectra to the resolving power(s) of the studied spectrum). 

3.2. List of reference lines 

The reference lines used to derive the surface gravity, iron abun- 
dance, andjndividual abundances were extracted from the line 
list of iBensbv et al.l d2003l) . We note that the astr ophysical oscil- 
lator strengths provided by IBensbv et al.1 d2003) were not used 
in calculating of the synthetic spectra (see Sect. [XT). For consis- 
tency, we used the Kurucz oscillator strengths for all lines (ref- 
erence or not). 



4. Performance assessment 

To assess the performances of SPADES, three series of tests 
were performed. 

The aim of the first series was to assess the internal system- 
atic and random errors obtained as a function of the S/N ratio, 
but independently of the question of the initialization parame- 
ters. We also investigated the correlations between the errors on 
the different parameters. 

In the second tests, high S/N ground-based observed spectra 
were analyzed with SPADES to assess the external errors. 

SPADES proceeds iteratively to derive the atmospheric pa- 
rameters. The third series of tests aims to assess the convergence 
properties of SPADES when it is initialized with parameters that 
are significantly off the true values. 

All tests were performed for the wavelength range and 
the resolving power of the HR13 (R = 22 500) and HR14B 
(R = 28 000) Giraffe setups. The three series of tests are pre- 
sented below. 



Result Teff distribution Result log(g) distribution 




5770 5790 5810 5830 4.2 4.3 4.4 4.5 4.6 4.7 4.8 

Teff (K) log(g) 



Result [Fe/H] distribution Result [Ni/Fe] distribution 




i i i i i i i i i i i 

-0.08 -0.04 0.00 0.02 0.04 -0.05 0.00 0.05 0.10 

[Fe/H] [Ni/Fe] 



Fig. 8. Distributions of the r eff (top left), log g (top right), [Fe/H] 
(bottom left) and [Ni/Fe] (bottom right), estimated by SPADES 
for 200 S/N=50 simulated HR13 and HR14B Giraffe spectra. 
The true parameters are T e g = 5800 K, logg = 4.5, [Fe/H] 
= 0.0 dex, and [Ni/Fe] = 0.0 dex 



4. 1 . Internal errors 

The internal systematic and random errors were assessed 
by Monte Carlo method. Starting from a synthetic spectrum 
from the reference grid with r e ff = 5800 K, logg = 4.5, 
[Fe/H] = 0.0 dex, £ = 1 kms" 1 and solar abundance scale, 
three series of 200 noisy spectra were generated for SNRs 
30, 50 and 100. The noise was injected following a Poisson 
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distribution. The 600 spectra were analyzed with SPADES 
to estimate their T e g, logg , [Fe/H] as well as [Ti/Fe] and 
[Ni/Fe] ratios as examples of abundance ratios estimates. In 
this first test, SPADES was initialized with the true parameters 
of the synthetic spectra. Table [2] presents the medians and 
dispersions of the distributions of residuals, i.e. estimated 
minus true, of the parameters derived by SPADES for each S/N 
ratio. Figure [8] shows the distributions of r e ff, logg , [Fe/H] 
, and [Ni/Fe] estimated by SPADES for the 200 S/N=50 spectra. 




4.0 4.2 4.4 4.6 4.8 5.0 

109(0) 



The distributions of temperature and gravity residuals 
present no significant biases. The systematic errors on the iron 
over hydrogen and titanium and nickel over iron ratios are of 
about a few milli-dex and are most of the time one order of mag- 
nitude smaller than the random errors. The random errors are 
small, e.g. at S/N=30: less than 30 K for T eS , 0.15 for log g and 
a few hundredths of dex for the abundance ratios, we emphasize 
that these are internal errors accounting only for photon noise. 
External errors are assessed in Sect. 14.21 



Table 2. Medians and dispersions of the distributions of resid- 
uals, i.e. estimated minus true, on T s g, logg , [Fe/H] , [Ti/Fe] 
and [Ni/Fe] derived by SPADES for the S/N: 100, 50 and 30. 
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Figure [9] presents the estimate by SPADES of T e g, log g , 
and [Fe/H] versus each other for the 200 stars with S/N=100 
(top), S/N=50 (middle), and S/N=30 (bottom). The clearer er- 
ror correlation is between the effective temperature and the iron 
abundance: an error of -50 K corresponds to an error of about 
-0.075 dex. This correlation, as the others, is intrinsic to the 
physics of the stell ar atmopsheres. While H„ is weakly s ensi- 
tive to metallicity dvan't Veer-Menneret & M egessier 1996j), the 
equivalent width of the iron lines is sensitive to the effective tem- 
perature. 



4.2. External errors: ground-based observed spectra 

The first series of tests was devoted to internal errors. The source 
of noise was Poissonian photon noise. The aim of the second 
series was to assess the external systematic errors. There are 
several possible sources of external systematics. In particular, 
the physics of the stellar atmospheres do not reproduce the 
reality perfectly (e.g. discrepant oscillator strengths, local 
thermodynamical Equilibrium assumption, ID modeling of the 
atmosphere), there are usually instrumental calibration residuals 
(e.g. from bias, flat-field, possibly fringing, etc.) and residuals 
from the observing conditions (e.g. from sky-subtraction). 
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Fig. 9. SPADES estimate of T e s, logg , and [Fe/H] versus each 
others for the 200 stars with S/N=100 (top), S/N=50 (middle), 
and S/N=30 (bottom). 



4.2.1 . The Giraffe solar spectrum 

To investi gate these problems, w e retrieved from the Giraffe 
archive (iRover et all I2005L I2008I) the high S/N ratio HR13 
(R = 22 500) and HR14B (R =28 000) observed solar spectrum. 
The Giraffe solar spectrum was acquired in twilight conditions , 
which could alter the equivalent widths (Mol aro et al.l [2008 ). 
Moreover, the library of reference spectra used by SPADES was 
calulated with the Kurucz oscillator strengths (log gf), not with 
astrophysical log gf, even for the reference lines used by the dif- 
ferent diagnostics. We therefore proceeded in two steps to ana- 
lyze the Giraffe solar spectrum. First the effective temperature 
and surface gravity were fixed to their bibliographic values, i.e. 
T e ff = 5777 K and log g - 4.44, and only the iron over hydro- 
gen and titanium and nickel over iron ratios were determined 
by SPADES. Then, in a second step, the solar spectrum was re- 
analyzed, this time requesting SPADES to determine the five pa- 
rameters r e ff, logg , [Fe/H] , [Ti/Fe] , and [Ni/Fe] . 

Table [3] lists the iron over hydrogen and titanium and nickel 
over iron ratios derived by SPADES when the effective tempera- 
ture and surface gravity were fixed to their bibliographic values 
Teg = 5777 K and log g = AAA. The iron over hydrogen and tita- 
nium over iron ratios are underestimated by -0.08 and -0.09 dex 
respectively, while the nickel over iron ratio is close to the true 
value. 

The solar spectrum was then re-analyzed, requesting 
SPADES to also determine T^ff and log g . SPADES was initial- 
ized with the bibliographic parameters Sun and converged in five 
iterations. Table |4]presents the parameters obtained by SPADES 
for the Giraffe solar spectrum. 

The effective temperature derived by SPADES, 5689 K, is 
88 K lower than the bibliographic temperature T e g = 5777 K. 
Figure [10] shows a zoom on the red wing of the H„ line 
from the Giraffe HR14B solar spectrum (black line). Two syn- 
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Table 3. iron over hydrogen and titanium and nickel over iron 
ratios derived by SPADES for a solar HR13 and HR14B so- 
lar spectrum, when the effective temperature and surface grav- 
ity were fixed to their bibliographic values r e jj = 5777 K and 
logg = 4.44. 



[Fe/H] 

(dex) 



[Ti/Fe] 

(dex) 



[Ni/Fe] 

(dex) 



-0.08 



-0.09 



-0.03 



Table 4. Effective temperature, surface gravity, iron over hydro- 
gen and titanium and nickel over iron ratios derived by SPADES 
for the Giraffe HR13 and HR14B solar spectra. 



Teff 

(K) 



[Fe/H] 
(dex) 



[Ti/Fe] 
(dex) 



[Ni/Fe] 
(dex) 



5689 4.46 



-0.18 



-0.06 



0.04 



thetic spectra are overplotted: in red a synthetic spectrum with 
r eff = 5689 K and in blue r eff = 5777 K (the values of the 
other parameters are those of the Sun: \ogg = 4.44 and [Fe/H] 
= 0.0 dex). The red spectrum (T e s = 5689 K) agrees much better 
agreement with the Giraffe solar spectrum than the blue spec- 
trum (r eff = 5777 K). 
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Fig. 10. Zoom on the red wing of the H a line from the Giraffe 
HR14B solar spectrum (black line). Two synthetic spectra are 
overplotted: in red a synthetic spectrum with T e fi = 5689 K and 
in blue T e ff = 5777 K (the values of the other parameters are 
those of the Sun: log g = 4.44 and [Fe/H] = 0.0 dex). 

The surface gravity derived by SPADES, logg = 4.46 is 
consistent with the Sun's gravity log g = 4.44. 

The iron to hydrogen ratio derived by SPADES is underes- 
timated by -0.18 dex with respect to the Sun. There are two 
main origins for this offset. On the one hand, there is the offset 



of -0.08 dex reported in table [3] that is likely due to the mod- 
eling of the spectra, in particular slightly discrepant oscillator 
strengths, and also to possible equival ent width alteration s due 
to twilight observing, as reported by Mol aro et alJ ([2008). On 
the other hand, as discussed in Sect. 14.11 the error on the effec- 
tive temperature propagates into an error on the iron abundance. 
The offset of -88 K here propagates into an offset of -0.1 dex. 

The titanium and nickel over iron ratios derived when 
r e (jand logg are determined by SPADES are similar to those 
obtained when r e ftand logg were fixed to their bibliographic 
values. The reason is that the iron, titanium and nickel abun- 
dances are all similarly affected by the -88 K offset and 
therefore the effect partly canceled out in the logarithmic ratios 
of titanium over iron and nickel over iron. 



4.2.2. S4N reference spectra 

The parameters of the Sun are very accurately known, which 
makes it a natural t arget to asses s the pe rformances of SPADES. 
On the other hand, Mol aro et al.l ((2008) reported an alteration of 
the equivalent widths of the lines in the Giraffe solar spectrum, 
which was acquired in twilight conditions. 

To further test the performances of SPADES with a ground- 
based observed spectrum, we selected 5 stars, the Sun, v And 
(HD9826), p Vi r (HP 102870) n H er (HD161797) and a Dra 
(HD 185 144) in ICavreletafl (1201 lh . The authors derived the 
effective temperatures of their stars from angular diameters 
measured with interferometry. All their stars have diameters 
known to 2% or better, leading to an accuracy on the effective 
temperatures o f about 25 K . All st ars are bright and were exten- 
sively studied. ICavrel et all (1201 ll) adopted for lo gg and [Fe/H] 
recen t determinations from the PASTEL catalog (Sou biran et al.l 
2010). We proceeded similarly, adopting for logg and [Fe/H] 
the median over the five most recent determinations contained 
in PASTEL. The titanium and nicke l over iron ratios were 
extracted from lValenti & Fischer! ([2005 ). Th e spectra for the five 
stars were retrieved from the S4N library (Allende Prieto et all 
2004). They were degraded in resolving power and in sampling 
and restricted in wavelength to the Giraffe HR13 and HR14B 
setups. 



SPADES converged in two (v And) to four (fi Her and <x Dra) 
iterations. Table |5]lists the bibliographic parameters and the pa- 
rameters derived by SPADES for the five tests stars. Fig.fTTIdis- 
plays the temperatures, gravities, iron over hydrogen and nickel 
over iron ratios estimated by SPADES versus their bibliographic 
values. Table |6]presents the mean and dispersion of the residuals 
(SPADES minus bibliographic) on the estimation of the atmo- 
spheric parameters as well as the titanium and nickel over iron 
ratios. 

On average, the effective tempe rature is recovere d with a 
bias of +46 K. For the same stars, ICavrel et al.l (1201 lh found a 
bias of similar amplitude, but with the opposite sign, of about 
-90 K. In addition to the small statistics and the details of 
the H ff profile fiting, we note that the physics used to com- 
pute the grids of models and synthetic spectra differ in several 
aspects. SPADES relies on Kurucz ATLAS9 models computed 
with the mixing length parameter 1/H ? = 1 .25 and SYNTHE 
spectra (see Sect. 13. U . ICavrel et alJ (1201 lh used the Kurucz 
ATLAS9 models with 1/H p = 0.5 and a modified version of 
the Kurucz BALMER9 with the Stark broadening treatment of 
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Table 5. Bibliographic parameters and parameters derived by 
SPADES for the five tests stars. 
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Fig. 11. Temperatures (top left), gravities (top right), iron over 
hydrogen (bottom left), and nickel over iron ratios (bottom right) 
estimated by SPADES versus their bibliographic values. 

Table 6. Means and dispersions of the residuals (SPADES minus 
bibliographic) on the estimation of the atmospheric parameters 
as well as the titanium and nickel over iron ratios. 



(K) 



[Fe/H] 
(dex) 



[Ti/Fe] 
(dex) 



[Ni/Fe] 
(dex) 



mean 

(T 



46 
25 



-0.08 
0.37 



0.01 
0.05 



-0.05 
0.12 



0.01 
0.09 



4.3. Convergence properties 
4.3.1. First test 

SPADES proceeds iteratively to derive the atmospheric parame- 
ters r e ff, logg , and [Fe/H] . In the first series of tests, SPADES 
was initialized with the true parameters. Of course, in reality, the 
true parameters are unknown. The aim of this third series of tests 
is to assess the behavior of SPADES when it is initialized with 
parameters that are significantly different from the true values. 

For this third test, 34 single stars were generated. Their pa- 
rameters were randomly drawn following uniform distributions 
over the following intervals: r e ffin [5200, 6400] K, logg in [3.0, 
4.5], and [Fe/H] in [-2.0, 0.6] dex. No star was simulated for 
[Fe/H] below -2.0, because the four Fe II lines present in the 
Giraffe HR13 and HR14B become too weak to allow derivat- 
ing the surface gravity. This is illustrated in the top right plot 
of Figure [12] which shows three synthetic spectra of metal- 
licity [Fe/H] = 0.0 dex (blue), [Fe/H] = -1.0 dex (red), 
[Fe/H] = -2.0 dex (black). The line at 6456.39 A is an Fe II 
line. At [Fe/H] = -2.0 dex, the line is barely visible. 
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Fig. 12. Extract of three synthetic spectra of metallicity [Fe/H] 
=0.0 dex (blue), [Fe/H] =-1.0 dex (red), [Fe/H] =-2.0 dex 
(black), showing the influence of the metallicity on an Fe I line 
(top left), an Fe II line at 6456.39 A (top right), a Ni I line (bot- 
tom left) and a Ti I line at 6258.1 1 A (bottom right). 



1999b and impact broadening by neutral hy- 



Alla rdet al.l (l2008i. 



IStehle & Hutcheonl 
drogen collisions of 

The iron over hydrogen and nickel over iron ratios are 
recovered without significant offsets. The titanium over iron 
ratio, determined from four Ti lines (without astrophysical 
log gf adjustement), presents an offset of -0.05 dex. To the limit 
of the small statistics, the error on r e ffShows no trend with T e g . 



For each of the 34 single stars, three noisy HR13 and HR14B 
Giraffe spectra were simulated for S/N of 30, 50 and 100. The 
noise was Poissonian photon noise. The three series of 34 spec- 
tra were analyzed by SPADES. In all 102 cases SPADES was 
initialized with the same parameters: T e ff = 5800 K, logg = 3.0 
and [Fe/H] = -1.0 dex. The left side of Figure [13] presents the 
differences between the true parameters of the 34 synthetic stars 
and the initialization parameters: r e ff(top), logg (middle) and 
[Fe/H] (bottom). 
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Table 7. Medians and dispersions of the distributions of T e g, 
\ogg , and [Fe/H] residuals, i.e. estimated minus true, derived 
by SPADES for the SNRs 100, 50 and 30. 
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Fig. 13. Left: Differences between the true parameters of the 34 
synthetic stars and the initialization parameters: r e ff(top), \ogg 
(middle) and [Fe/H] (bottom). Right: Distributions of r e ff(top), 
logg (middle) and [Fe/H] (bottom) residuals for S/N=50. 



SPADES converged for 33 out of the 34 stars in at most 
six iterations (for all SNRs). It did not converge for one star 
with the parameters: r e ff = 5386 K, \ogg — 4.3, and [Fe/H] 
= 0.4 dex. At the first iteration, the difference of 500 K between 
the true r e ff and the initialization value and the high true metallic - 
ity [Fe/H] = 0.4 dex, led SPADES to look for a synthetic spec- 
trum with a metallicity higher than 1 dex, i.e. beyond the upper 
boundary of the library of synthetic spectra. Unable to find the 
synthetic spectrum, SPADES stopped processing this spectrum 
(at all three SNRs). 

In all 33 converging stars, it was observed that the effective 
temperature converged faster than the other parameters. This is 
because in mid-F-G dwarfs the H a wings are hardly sensitive 
to wr ong \ogg and [Fe/H] (Ivan't Veer-Menneret & Megessierl 
1 19961) . while on the contrary \ogg and [Fe/H] are sensitive 
to wrong r e ff. This suggests a possible future evolution for 
SPADES, which could at the first iteration determine only 
r e ff and then at the next iterations determine all three atmospheric 
parameters. This would probably have prevented SPADES from 
looking for spectra outside the boundaries of the library for the 
34'' 1 spectrum. 

Table |7]presents the medians and the dispersions of the dis- 
tributions of r e ff, logg , and [Fe/H] residuals, i.e. estimated mi- 
nus true, derived by SPADES for each S/N ratio. The right side 
of Figure Qj] shows the distributions of r e ff(top), logg (middle), 
and [Fe/H] (bottom) residuals for S/N=50. The performances 
derived from test 3 (for of a wider variety of input atmospheric 
parameters) are similar to those obtained in test 1 . The conver- 
gence process has therefore a very moderate impact on SPADES 
final r e ff, \ogg , and [Fe/H] . 

Out of the 33 converged stars, [Ti/Fe] and [Ni/Fe] ratios 
were derived for 2 1 . Below ametallicity of about - 1 .4 dex (the 
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precise value is a function of r e ff), the lines contained in the 
HR13 and HR14B become too weak to derive the abundances. 
This is illustrated by the bottom left and bottom right plots in 
Figure[T2] At [Fe/H] — -1 dex (red line), the Ni and Ti lines (at 
6258.11 A) are already weak and at [Fe/H] = -2 dex they are 
not visible anymore. 

Table [8]presents the medians and the dispersions of the dis- 
tributions of [Ti/Fe] and [Ni/Fe] residuals, i.e. estimated minus 
true, derived by SPADES for the SNRs 100, 50 and 30. The sys- 
tematic and random errors derived from test 3 are on average 
larger than those obtained from test 1 . The main reason is very 
likely that in test 3, [Ti/Fe] and [Ni/Fe] ratios are obtained for 
a wide range of metallicities down to [Fe/H] ~ -1.4 dex, while 
in test 1 a solar abundance star was considered. 



Table 8. Medians and dispersions of the distributions of [Ti/Fe] 
and [Ni/Fe] residuals, i.e. estimated minus true, derived by 
SPADES for the SNRs 100, 50 and 30. 
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4.3.2. Second test 

Iterative methods present the risk to converge toward a local 
minimum, rather than the global minimum. To investigate this 
risk in more detail, we selected one out of the 34 times 200 spec- 
tra analyzed in Sect. 14. 3. TI above, with S/N=50. The true param- 
eters of the spectrum were T e $ = 5992 K, logg = 4.4, [Fe/H] 
= -0.3 dex, [Ti/Fe] = dex and [Ni/Fe] = dex. SPADES an- 
alyzed the spectrum a first time, starting from its true parameters 
and found T eff = 6025 K, logg = 4.35, [Fe/H] = -0.26 dex, 
[Ti/Fe] = 0.04 dex and [Ni/Fe] = 0.05 dex. SPADES then re- 
analyzed the spectrum another 50 times, but with starting pa- 
rameters randomly chosen following uniform distributions cen- 
tered on the true parameters and widths of + 500 K (r e ff), + 0.5 
(logg ) and + 0.4 dex ([Fe/H] ). The volume in the parameters 
space considered for the initialization is larger, i.e. conservative, 
with respect to the estimated performances of SPADES-TE (see 
Sect. 12.5.2b . Table |9]presents the means and residuals of the pa- 
rameters determined by SPADES when initialized randomly mi- 
nus those determined by SPADES when initialized with the true 
parameters. Means and dispersions are all low, showing that at 
S/N=50, SPADES is very weakly sensitive to the initialization 
conditions. 
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Table 9. Means and dispersions of the residuals (SPADES ini- 
tialized randomly minus SPADES initialized with the true pa- 
rameters) on the estimation of the atmospheric parameters as 
well as the titanium and nickel over iron ratios. 

T cff log# [Fe/H] [Ti/Fe] [Ni/Fe] 

(K) (dex) (dex) (dex) 

mean 1 -0.01 -0.007 0.002 -0.007 

cr 3 0.02 0.007 0.004 0.007 



5. Summary 

The automated stellar parameters determination software 
SPADES was presented. It relies on a line by line comparison 
between the studied spectrum and synthetic spectra. SPADES 
derives the radial velocity, the effective temperature, the surface 
gravity, and the iron and individual abundances. 

SPADES internal systematic and random errors were 
assessed by Monte Carlo method. For example, simulating 
Giraffe HR13 and HR14B setups from synthetic spectra for 
a star with r eff = 5800 K, logg = 4.5, [Fe/H] = 0.0 dex 
and for a S/N of 100, the stellar parameters were recovered 
with no significant bias and with 1-cr precisions of 8 K for the 
temperature, 0.05 for the logg , 0.009 for [Fe/H] , 0.003 for 
[Ti/Fe] , and 0.01 for [Ni/Fe] . The external systematic errors 
were evaluated with ground-based observed spectra and the 
convergence properties were assessed by a Monte Carlo method. 

Several evolutions of SPADES are planned. On the one hand, 
we aim to adapt and test SPADES for a greater variety of spec- 
tral types and luminosity classes. This includes addingare the 
micro-turbulence to the parameters determined by SPADES. On 
the other hand, we plant to implement additional diagnostics to 
determine the atmospheric parameters, e.g. use the wings of the 
strong lines to derive the surface gravity (for example Mg lb 
triplet in dwarfs or infra-red Ca II triplet in giants). 
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